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We extend the "analogue spacetime" programme by investigating a condensed-matter system that 
is in principle capable of simulating the massive Klein-Gordon equation in curved spacetime. Since 
many elementary particles have mass, this is an essential step in building realistic analogue models, 
and a first step towards simulating quantum gravity phenomenology. Specifically, we consider the 
class of two-component BECs subject to laser-induced transitions between the components. This 
system exhibits a complicated spectrum of normal mode excitations, which can be viewed as two 
interacting phonon modes that exhibit the phenomenon of "refringence" . We study the conditions 
required to make these two phonon modes decouple. Once decoupled, the two distinct phonons 
generically couple to distinct effective spacetimes, representing a bi-metric model, with one of the 
modes acquiring a mass. In the eikonal limit the massive mode exhibits the dispersion relation of 
a massive relativistic particle u> — + c 2 k 2 , plus curved-space modifications. Furthermore, it is 
possible to tune the system so that both modes can be to arranged travel at the same speed, in which 
case the two phonon excitations couple to the same effective metric. From the analogue spacetime 
perspective this situation corresponds to the Einstein equivalence principle being satisfied. 
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I. INTRODUCTION 

Analogue models for gravitation (which should more 
accurately be referred to as analogue models for curved 
space-time) can be used to simulate classical and quan- 
tum field theory in curved space-time [H, IE IE 0, IE 

IlEIIEEIIlEEEElSEEEEliEliJ The 

first analogue model for black holes, and for simulating 
Hawking evaporation, was suggested by Bill Unruh [ij in 
1981. He demonstrated that a sound wave propagating 
though a converging fluid flow exhibits the same kine- 
matics as does light in the presence of a curved space- 
time background. Since then several other media - 
e.g. flowing dielectrics and quantum liquids Q - 
have been analyzed, and the field has developed tremen- 
dously EE El E3- 

The first approach specifically us- 
ing Bose-Einstein condensates as an analogue model was 
made some nineteen years after Unruh's original pa- 
per . Since then various configurations of BECs have 
been studied to simulate different scenarios for grav- 
ity |E IE IE IE ■ Until now it has only been possible 
to simulate photons, (generally speaking, massless rela- 
tivistic particles) , propagati ng through a curved space- 
time El EE EE El El HE El El E3 In the present 
article a two-species BEC is used to extend the class of 
equations that can be simulated to the full curved-space 
massive Klein-Gordon equation in (3+1) dimensions. (In 
the language of the BEC community this corresponds to 
a specific and technologically interesting way of giving a 
mass to the phonon.) From the viewpoint of the gen- 
eral relativity community, this article provides a route 
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for analogue simulations of curved-space quantum field 
theory that are more general and realistic than those 
considered to date. (A preliminary sketch of these of 
these ideas appeared in reference [l9j. An alternative 
route to the Klein-Gordon equation via textures in 3He 
is described in |2(|, and the use of "wave guides" to re- 
strict transverse oscillations and provide an effective 1+1 
dimensional mass is described in 9J.) 

II. TWO-COMPONENT BECS 

The class of system we will use in our theoretical anal- 
ysis is a two-component BEC. More specifically we con- 
sider an ultra-cold two-component BEC atomic gas such 
as, for example, a two-component condensate of 87 Rb 
atoms in different hyperfine levels, which we label \A) 
and \B). (Experiments using two different spin states, 
\F = 1, m = —1) and \F = 2,m = 2), were first per- 
formed at JILA in 1999 21].) At very low temperatures 
nearly all atoms occupy the ground state. For the follow- 
ing calculation the non-condensed atoms are neglected. 
(In HE EE El EE El finite temperature effects are taken 
into account.) The quantized field describing the micro- 
scopic system can be replaced by a classical mean-field, 
a macroscopic wave-function. In this so-called mean- 
field approximation the number of non-condensed atoms 
is small. Interactions between the condensed and non- 
condensed atoms are neglected in the mathematical de- 
scription, but two-particle collisions between condensed 
atoms are included. In the case of a two-component sys- 
tem, interactions within each species {Uaa, Ubb) and 
between the different species {Uab = Uba) take place. 
In addition the two condensates are coupled by a laser- 
beam, which drives transitions between the two hyperfine 
states with a constant rate A. Without the laser coupling 
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A, no mass term is generated, which is consistent with the 
analysis in reference 0. (In that article the advantages 
of using a two-component BEC to simulate cosmological 
inflation are presented.) Because, for current purposes, 
the two species are different hyperfine levels of the same 
atom, the masses of the individual atoms are approxi- 
mately equal (to about one part in 10 16 ). In the current 
article we therefore set 



rriA = m B 



m, 



(1) 



though we note that strictly speaking ttla ^ tub- (Indeed 
in more abstract situations, potentially of relevance in 
quantum gravity phenomenology, it is worthwhile and 
instructive to keep this extra mass dependence explicit.) 

The resulting coupled two-component time-dependent 
Gross-Pitaevskii equations are: 



— V 2 +V A - HA 
2m 



-U AA \^a\ 2 + U AB \^b\ 2 



(2) 



ihdt^s 



2 m 



V + Vb ~ Pb 



(3) 



-U BB \^b\ 2 + U A b\^a\ 2 



1 $B + W A 



where Va/b denotes the two extern al p otentials, and 
Pa/b the t w0 chemical potentials [2?! |28|. We note that 
the parameter A can be either positive or negative with- 
out restriction, while the interaction parameters Uaa, 
Uab, and Ubb are typically though not always positive. 
Adopting the Madelung representation, the two conden- 
sate wave- functions ^ x can be described in terms of their 
densities {pa,Pb} and phases {9a, 9b}- 



x 



Ipx e 



for X = A, B . 



(4) 



These four variables in general depend on both time and 
space. 



III. WAVE EQUATION 



where po(r) is the density of the ground state, a is the 
scattering length, m is the atomic mass and U is the self- 
interaction constant. The mathematical equations de- 
scribing these perturbation lead to the well-known hydro- 
dynamic equations, which are the basis for the most fruit- 
ful of the analogies between condensed matter physics 
and general relativity P, 0, flfl ITfij . 

An extension of this method can be used to obtain the 
kinematic equations for small perturbations propagating 
in a two-component system. Given that the density mod- 
ulation is small, the perturbations in the densities and 
phases can be expanded around their macroscopic states 
ijljl using perturbation theory: 



= V 'Pxo +ep X ie 



i(6 xo +s8 xl ) 



for X — A,B . (6) 



These states still satisfy the coupled Gross-Pitaevskii 
equation. When developing a perturbation analysis for 
the fluctuations we find that unless we set the back- 
ground phases equal to each other (9aq = 6 bo) the 
calculation becomes quite intractable. Specifically, one 
encounters damping and damping terms dependent on 
Aab = 9 ao — 9 bo- In the appendix (|A"I) we briefly present 
the result obtained for arbitrary — even time-dependent 
- background phases. While these terms and their im- 
plications are of interest in their own right, in the follow- 
ing focus will be set on Aab = automatically implying, 
in particular, that the background velocities of the con- 
densates, 

vao = {h/m)V9 A o and v B o = {ft/m)V6 B Q, (7) 
are equal: 



Vo 



VAO = VBO- 



(8) 



After a straightforward calculation, the terms of first 
order in e include two coupled equations for the pertur- 
bation of the phases 



a - Y7A Uaa Uab 

7 AX = - v ■ v9 A i z— Pax r~ Pbi, 

n n 

n ™ Ubb Uab 
TB\ =- v ■ Vf/Bi r— Pbi r— Pax- 



(9) 



We study zero sound in the overlap region of the two- 
component system, produced by exciting density pertur- 
bations which are small compared to the density of each 
condensate cloud. In the first experiment studying local- 
ized excitations in a one-component BEC j2£j , the optical 
dipole force of a focused laser beam was used to modify 
the trapping potential, generating a small density mod- 
ulation. Using phase-contrast imaging it was shown that 
the resulting perturbation corresponds to a sound wave. 
The observed speed of sound is 



Here 



c(r) 



47T h 2 a po{f) 



Upo(r) 



(5) 



Uaa 
Ubb 
Uab 



= Uaa-± 



Ubb 



Uab 



2 {pAof 2 ' 

A y/PAO 
^(PBO) 3 / 2 ' 

A 1 

2 y/PAO Pbo ' 



(10) 



are modified interaction potentials for the two coupled 
condensates. In addition to these two phase equations, 
there are two coupled equations for the density pertur- 
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bations 



Pai 



PB1 



- V pao^Oai + PAI Vq 

2A 

+ -^Vpao Pbo {Obi - &ai), 

- V f — pbq^Obi + PBl v 
2A 



(11) 



+ -r- \fp~M) Pbo (9ai ~ Obi)- 

To adopt a more compact representation of the physics, 
it is useful to define several matrices and vectors. First, 
define the coupling matrix 



Uaa Uab 
Uab Ubb 



A second coupling matrix, defined as 



A = 



2A y/pAO PBO 



+1 -1 
-1 +1 



(12) 



(13) 



vanishes completely if the laser coupling A is switched 
off. Last but not least, it is also useful to introduce the 
mass-density matrix D 



D = 



PAQ 
PBO 



Now define the two-component column vectors 



0A1 
Obi 



and 



P = 



PAY 

Pbi 



(14) 



(15) 



Collecting terms into a 2 x 2 matrix equation, the equa- 
tions for the phases ijUJl and densities l|ll(l become a com- 
pact pair of first-order PDEs 



= - S p - v ■ V9, 



-V ■ (£> V0- 



pvo 



A0. 



(16) 



(17) 



Equation i|16f) can be used to eliminate p and p in 
equation H17[) . leaving us with a single second-order 2x2 
matrix equation for the perturbed phases: 

d t (s- 1 0) = -d t (E' 1 v - V0) - V- (v S- 1 0) 

+V- [(D-v E- 1 v -) V0] +A0. (18) 

This equation tells us how a localized collective excitation 
in a A-coupled two-component BEC develops in time. It 
is a special case of the "normal mode" formalism devel- 
oped in [n|,H0. 

If we adopt a n-dimensional "spacetime" notation by 
writing x a = (£,£*), with i S {1, 2, ••• ,n — 1} and 



a e {0, 1, 2, • • • , n — l}, then this equation can be very 
compactly rewritten as 0, : 



d a (f ab d b 9) + A = 0. 



(19) 



Here f ab is a n-dimensional rank 2 tensor density, each 
of whose components is a 2 x 2 matrix. Specifically 



00 _ rs—l, fOi _ n — 1 „,i _ riO 



r = — h vi = r 



and 



DSi. 



(20) 



(21) 



So far there are no significant restrictions on the back- 
ground densities {pao, Pbo), the joint background flow 
velocity vq, the interaction constants (Uaa, Ubb, Uab), 
and coupling constant, A. If we do not decouple the two 
modes, then the most we can say is that the coupled 
system exhibits "refringence" in the sense of [Hill. 



IV. MODE DECOUPLING 

The first step in analyzing equation l|18f) . or the equiv- 
alent (|19|) . is to ask whether it is possible to tune the 
system so as to completely decouple it into two inde- 
pendent phonon modes. Only if the two modes decou- 
ple is it possible to assign individual masses and space- 
time geometries to the decoupled modes 0, 0, 131. In 
the absence of decoupling, one simply has a complicated 
two-component system with no clear spacetime interpre- 
tation. In performing the analysis we have found that 
decoupling is not possible without introducing several 
constraints on the background quantities. 

The decoupling analysis can be performed in several 
different ways, all ultimately leading to qualitatively sim- 
ilar physics, with minor technical differences. The ma- 
jor decision to be made is whether one imposes decou- 
pling at the level of physical acoustics (at the level of 
the wave equation) or at the level of geometrical acous- 
tics (at the level of dispersion relations). The fact that 
physical acoustics leads to propagation phenomena more 
subtle than those detectable in the geometric acoustics 
limit is well known 0, 0] . A particularly illustrative 
example is provided by acoustic propagation in a fluid 
with non-zero vorticity [30| . where the geometric acous- 
tics approximation leads directly to a conformal class of 
effective spacetime metrics, while the physical acoustic 
approximation leads to a complicated system of PDEs. 
If the only thing you can detect experimentally is the 
dispersion relation, then one should adopt geometrical 
acoustics and not demand to decouple the wave equa- 
tion itself. On the other hand, if one has experimental 
probes that couple directly to the wave itself, then the 
decoupling should be performed at the level of physical 
acoustics. We shall do both, and compare results later. 
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V. PHYSICAL ACOUSTICS 

At the level of physical acoustics one treats the wave 
equation (|19|) as primary, then decoupling requires that 
all the (symmetric) matrices f ab , and the (symmetric) 
matrix A, be simultaneously diagonalizable by position- 
independent orthogonal matrices O. That is, decoupling 
requires 



this condition fails and S is singular. Then equation H6|) 
cannot be solved for the column vector p and we cannot 
even set up the wave equations (|T%|) or (|T§|| . ] All in all 

we must have 



Uaa Uab 
Uab Uaa 



det(S) ^ 0. 



(29) 



jab = Q T f a^ Q . A = T A di: 

since then equation (|19|l becomes 



<:), 



'a ^/di 



Now since 



diag 



A 



d b e) +A diag e = o 



Secondly, we are now left with the penultimate term 
; Q — O 9; (22) m equation l (T%|l . Because the eigenvectors of S" 1 are 
already known, there is only the mass-density matrix D 
to consider. The final constraint to decouple the equation 
for the two phase perturbations is now easily seen to be 
(23) pao = Pbo- We shall simply denote this common density 
by pq. That is 



(24) 



the matrix that diagonalizes it is clearly position- 
independent, and the condition for simultaneous diag- 
onalization reduces to 



[f a \ f cd 



0: 



[/ j A] = 0; 



(25) 



where the commutators are to be interpreted in the sense 
of 2 x 2 matrix multiplication. That is, the matrices S, 
D, and A must all be simultaneously diagonalizable. 

We could now proceed by direct calculation of the three 
commutators 



[A A]; [A, 3]. 



(26) 



A perhaps more direct analysis can be performed directly 
in terms of equation l|18|) . Focusing on the last term 
in equation l|18|) . the eigenvectors for non-zero coupling 
A ^ are given by { [1, 1] T , [—1, 1] T }. The correspond- 
ing eigenvalues are {0, 4 X^paqPbq/K\- These eigenvec- 
tors (though not the eigenvalues) are fixed, position- 
independent, and indeed independent of any of the other 
physical variables. As a result the only way to decouple 
equation l(18jl into two independent phonon modes, in a 
position-independent manner, is to demand: 



(27) 



We now analyze equation l|18l) term by term with respect 
to this particular decomposition. 

Firstly, the term on the LHS, and the first two terms 
on the RHS in equation l|18fl . have the same eigenvec- 
tors as equation l|27|) if and only if Uaa = Ubb- The 
eigenvalues of S~ x corresponding to the eigenvectors 
{[1,1] T , [-1,1] T } are then 



" 1 " 




' -i ' 


1 


+ e 2 


i— 1 



(Uaa + UabY (Uaa -Uab) 



(28) 



This places another constraint on the interaction vari- 
ables: Uaa ±Uab- [Indeed, note what happens if 



D = dl=^-I. 



(30) 



In view of this last constraint the first coupling ma- 
trix simplifies considerably (and this was certainly not 
obvious when we started the analysis). We now have 



Uaa 
Ubb 
Uab 



= u 



BB 



A 

W 

Uab + =— , 
2po 



(31) 



and the first mode-decoupling constraint reduces to 
Uaa = Ubb- In order for the eigenvalues of to 
be well defined we need both Uaa + Uab ^ and 
A 7^ Po{Uaa — Uab), which is a mild easily satisfiable 
constraint. 



VI. BI-METRICITY 

Applying all this to equation l|18() one now sees how 
mode decoupling is equivalent to bi-metricity: One ob- 
tains two [dimension-independent] decoupled equations 
for the phonon modes described by the eigenstates of 
equation 127(1 : 



d a (f1 b W 



4Ap 



6 2I 6 Il for 7=1,2, (32) 



where 62 1 is the usual Kronecker delta. Here 



ft 



a b 



d 


\ — '■ 


i 




-v J 


cj 8^ — Vq v j 



(33) 



where the possibly distinct propagation speeds cr are de- 
fined in terms of the eigenvalues Sj of the matrix 2 by 



" 2 r =Sid 



djUAA + i-lYUAs) 
h 



(34) 
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that is 



2 _ Po{Uaa + {-IYUab) 

c i — 

m 



(35) 



It is important to note that decoupling by itself does not 
force the two propagation speeds to be equal — decou- 
pling in this context generically produces a two-metric 
model [bi-metricity] , and the demand that every phonon 
couple to a single effective metric [mono-metricity] is 
an additional independent condition. (In general the 
phonons arising from a system of N interacting BECs 
would, if they decoupled in the above manner, lead to an 
iV-metric model.) 

Introducing the (dimension-independent) natural os- 
cillation frequency 



2 _ 4A po c\ 



we can write 



d a (ff d b e 



u 1 S 2I 



4A m c\ 



for 7=1,2, 



(36) 



(37) 



where ujq now has the physical interpretation that it 
is the frequency at which a position-independent (zero- 
momentum) configuration oscillates. 

When converting the contravariant tensor densities ff 
into covariant spacetime metrics g ab one encounters a 
number of dimension-dependent factors |lfjl Il8j . such 
that naive application of the formalism leads to peculiar 
dimensionalities for physical quantities. The best way we 
have found for keeping track of these factors is to intro- 
duce space and time independent reference constants c* 
and cL and write 



lib 



f 



ah 



ci d t 



1 

ci 



-1 



so that 
d,(fT 9j 



(38) 



c* d uiq S 2 i z 



!, for 7=1,2. (39) 



C2 c 2 

Here c* and d* are any convenient but fixed reference 
values. They might be (for instance) the spatial average 
values of c 2 and d over the entire condensate, or they 
might be chosen in terms of the speed of light and other 
fundamental constants. They are introduced for conve- 
nience, and do not ultimately affect the physics we are 
discussing. 

Now introducing a pair of effective "spacetime metrics" 
by the identifications 

]f and 5/ = _!-_, (40) 



-9i dt 



det[<?? 



we can recast these wave equations equations as a pair 
of curved-space Klein-Gordon [massive d'Alembertian] 
equations 



™2 Jl 
phonon * 



(41) 



where as we shall soon see all quantities carry their stan- 
dard dimensionality 

After a brief algebraic calculation we find 



and so 



at 



c* d 
ci d* 



-ai = ci 



-n/(n-2) 



d 



ci d* 



n/(n-2) 



(42) 



-1 



(43) 

These quantities depend on the space-time dimension 
n 0, 0| in such a manner that 



aib 



c* d 



n/(n-2) 



Si d* 

Finally the mass-term is 





7 

-v J 


7 


5» 



m 



phonon 



C2 d t 



-2/{n-2) 



(44) 



(45) 



We note that this "mass" term can in principle depend 
on position, so some authors might prefer to refer to it 
as a "potential" term. There is no universal agreement 
on terminology regarding this point, but the bulk of the 
community would be happy to refer to this as a "position 
dependent mass" . Furthermore, if one focuses on the core 
of the BEC cloud, where all gradients are by definition 
zero (or small), then this mass term is guaranteed to be 
approximately constant. 

Now the two metrics g^ b , the inverse metrics gt, and 
the phonon mass mp honon , all depend on the normaliza- 
tion constants c* and g?*. This is as it should be, since 
c* and d* were introduced to give canonical dimensions 
to these quantities. In particular, since c* and <i* con- 
tribute to the overall conformal factor in the analogue 
spacetime metric, they set the overall scale with respect 
to which analogue "masses" are to be measured. How- 
ever, correctly formulated physical questions will depend 
only on parameters such as ujq, cj, and vq which are inde- 
pendent of these normalization constants. For instance, 
in the eikonal limit the dispersion relation for these two 
decoupled phonon modes is simply 



v Q -k) -c]k 2 =uj 2 6 2 i . 



(46) 



(A similar calculation, but restricted to a one-condensate 
system, where all variables are likewise allowed to be time 
and space dependent, but no mass term is present, has 
been presented in Q.) 

For the phase vector 9±, (corresponding to perturba- 
tions in the two condensates A and B oscillating "in 
phase"), the mass term is always zero. However, for a 
laser-coupled system (A ^ 0) the mass-term in the equa- 
tion for 9 2 , (corresponding to perturbations in the two 
condensates A and B oscillating in "anti-phase"), does 
not vanish. 
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VII. MONO-METRICITY 

Comparing the definition for the speed of sound (J5J 
in a one component system, with the two speeds cj in- 
troduced here, we see that the cj are simply the 
A-modified speeds of sound for each decoupled phonon 
mode. (For instance, consider an idealized situation in 
which the two condensates decouple completely, Uab = 
and A = 0, the two c/'s become the independent phonon 
speeds in each condensate cloud.) This fact leaves us 
with the possibility of constructing two different types of 
analogue model. So far we have been dealing with a two- 
metric structure, which is interesting in itself [l2l l3l| . 
For instance, in the absence of laser-coupling, A = 0, the 
presence of two different speed of sounds can be used for 
tuning effects 0. 

If we wish to more accurately simulate the curved 
spacetime of our own universe, another constraint should 
be placed on the system, to make the two speeds of sound 
equal c — ci — C2 . This yields a single sound-cone struc- 
ture, to match the observed fact that our universe ex- 
hibits a single light-cone structure. This condition is ful- 
filled if we set 



Uab = 0; 
In this case 



A 



-2 Uab Po- 



Po U A a _ Po (U A a + U A b) 



(47) 



(48) 



while we have the dimension-independent result that the 
natural oscillation frequency becomes 



Po U A b(Uaa + Uab) 



h 2 



(49) 



In counterpoint 



8 Uab Po m 



phonon 



(u. 



AA 



U A i 



h 2 



ci po 



l/(i- 2 ) 



(50) 

We note that in typical situations U aa ~ Ubb and 
Uab are both positive, corresponding to repulsive atomic 
interactions |2g. This implies that A is then negative, 
corresponding to a negative trapping potential, but a 
positive ujq and a positive mJ hono]1 . 

It is also possible to choose systems such that Uab 
is negative, as long as c 2 is still kept positive, which in 
turn requires Uaa + Uab to be positive, and places a re- 
striction on the relative sizes of the self-interactions and 
cross-interactions. This implies that A is positive. Fur- 
thermore the natural oscillation frequency cu 2 and the 
phonon mass m p honon 2 will then be real and negative. 
That is, attractive atomic interactions, which signal an 
instability in the condensate, would in our analysis cor- 
respond to a tachyonic phonon. 

If in contrast we permit c 2 to go negative (that is 
we permit Uaa + Uab to be negative) then we have 
a "Euclidean signature regime" where phonons do not 



propagate — this corresponds to the gross instability of 
the condensate which manifests itself as the "Bose-Nova" 
phenomenon |32j |. 

In this mono-metric situation, while the in-phase per- 
turbations will propagate exactly at the speed of sound, 



vo + k c, 



(51) 



the anti-phase perturbations will move with a lower 
group velocity given by: 



duj 
dk 



—~ = vq + k 



i/Wq + ° 2 fc2 



(52) 



Here k is the usual wave number. This explicitly demon- 
strates that the group velocity of the anti-phase eigen- 
state depends on the laser-induced coupling between the 
condensates. 



VIII. SPECIAL CASE: CONSTANT S 

There is a special case that is worth considering sepa- 
rately. Suppose that the matrix 5 is a time and space- 
independent constant. Then by defining 



3 - 
'new i — 1 



(53) 



multiplying equation 1)16(1 by S +1 / 2 , and appropriately 
commuting 5 through the space and time derivatives, 
one can re- write equation l|16H as 

<9 2 £new = ~d t (I Vq ■ V0 now ) - V • [VQ I 9 nC w J 

+V ■ [{C 2 - v I v • ) W ncw ] - M 2 9 ncWl (54) 

where 

C 2 = ~}' 2 DE 1 / 2 ; M 2 = -~}' 2 k~}' 2 . (55) 

Now the existence of the matrix square root S 1 / 2 fol- 
lows since 3 itself is real and symmetric. Indeed by the 
Hamilton-Cayley theorem for 2x2 matrices we know 



S 1 / 2 = aI + &S, 



(56) 



where a and b are functions of the eigenvalues of S. A 
little matrix algebra yields (for 2x2 matrices only) the 
explicit formula 



ctV2 _ 



VdctS 1 + 5 



2VdetS + tr(S) 



(57) 



It is now clear that both C 2 and M 2 are symmetric ma- 
trices, so in this constant-2 special case mode decoupling 
of the wave equation l)54[l. that is bi-metricity, requires 
the simple constraint 



[C 2 , M 2 \ = 0, 



(58) 
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which is equivalent to the matrix equation 

DE A = AS D. (59) 

In terms of the physical parameters we obtain the con- 
straint 



Pao Uaa — Pbo Ubb = (pao — Pbo) Uai 



(60) 



Note that to get to this stage we had to assume S was 
constant (so Uaa, Uab, and Ubb are constant), and 
then from the above we see that decoupling requires that 
Pao / Pbo must be a position and time- independent con- 
stant. We do not however need pao = Pbo', by restricting 
the position and time dependence of S we have permit- 
ted other parts of the wave-equation to possess an al- 
gebraically more general solution to the decoupling con- 
straint (the bi-metricity constraint). Subject only to the 
condition ijfiUI) the speeds of sound and masses are sim- 
ply extracted as eigenvalues of the matrices C 2 and M 2 
respectively. 

Imposing the bi-metricity condition (|6U[1 the two eigen- 
values of C 2 , the two speeds of sound, are 



JLf 

2m I 



PaqUaa + PboUbb ± (pao + Pbo 



)Uab] 



(61) 

Additionally imposing the mono-metricity condition re- 
quires the two eigenvalues of C 2 to be the same. That 
is, mono-mericity enforces Uab — 0. Together with the 
decoupling condition equation (|60() . we now get: 



Pao Uaa = Pbo Ubb- 



That is 



2 Pao Uaa 
c = 



Pbo Ubb 



m 



m 



(62) 



(63) 



and 

ujI = tr[M 2 ] 

4 pao Pbo Uab 



Uaa + Ubb + Uab fe + ^f 



h 2 



(64) 



The eigenvector corresponding to the massless phonon 
mode is 9 ncw oc IvUbb, vUaa\ T , while that for the 



massive phonon mode is 6* now cx [—yuAA, \/Ubb] T ■ 
In terms of the original variables 9 this corresponds 
to decoupled eigenmodes 8 oc [1, 1] T for the mass- 
less mode, while that for the massive phonon mode 
9 oc [-Uaa, U B b] t ■ 



IX. SPECIAL CASE: CONSTANT D 

Similarly, consider the case where the matrix D is in- 
dependent of position and time — this corresponds to 



a situation where the background densities pao and pbo 
are constant, though they do not need to be equal. In 
this situation we can define 



9 - D +1/2 

"new ~~ ' J 



(65) 



Then multiplying equation (|16|l by D~ x / 2 , and appropri- 
ately commuting D through the space and time deriva- 
tives, one can re- write equation (|16|l as 

d t {C- 2 d t 9 ncw ) = -d t (C- 2 vo ■ V0 ncw ) 

-V ■ (Vo C- 2 9 ncw ) 

+V • [(I - v C- 2 v ■ ) W new ] - M 2 9 ncw , (66) 

where now 

c -2 = £-1/2 s _i ^-1/2. m 2 - -D-Va A D- 1 ' 2 . 

(67) 

The analysis is now similar to the case of constant-^. 
Decoupling (bi-metricity) requires 



that is 



[C- 2 ,M 2 ] = 0, 



D- 1 A = AD- 1 5" 



(68) 



(69) 



It is easy to rearrange this equation to yield 

A S D = D 3 A, (70) 

which is the same matrix equation as encountered in 
the constant-S case. Consequently, decoupling requires 
the same algebraic condition, equation (|60|) as in the 
constant-S case. The speeds of sound are now given by 
the eigenvalues of the matrix 



C 2 



£+1/2 s £+1/2 



(71) 



and explicit computation again yields, after imposing 
the decoupling constraint (|60|l . the same algebraic result 
(|61|l . Imposing mono-metricity again leads to Uab = 
and equation Q62f) . which in this case can be read off 
by inspection from equation l|71l) . In short, there is no 
new physics hiding in the case of constant D; we again 
see that by restricting the position and time dependence 
of some of the coefficients in the wave equation <|18[1 we 
have permitted other parts of the wave-equation to pos- 
sess a slightly more general solution to the decoupling- 
constraint (the bi-metricity constraint). 



X. GEOMETRICAL ACOUSTICS — FRESNEL 
EQUATION 

In geometrical acoustics we adopt the eikonal approx- 
imation 



A exp(-itp) 



(72) 
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with a slowly varying amplitude A and a rapidly varying 
phase ip. We also assume that the coefficients in the 
differential equation (|18JI are slowly varying compared to 
the phase. This approximation leads, before we apply the 
decoupling constraints, to the Fresnel equation [IlLll2| 



det {lo 2 "T 1 - 2uj (v Q ■ k) -T 1 

-[D k 2 -(v -k) 2 S- 1 ]+A]}=0. 



(73) 



In a 2-component BEC the Fresnel equation is in general 
a quartic dispersion relation for two interacting phonon 
modes. This approximation makes sense when the period 
and wavelength of the phonon mode is small compared 
to the timescale and distance scale over which the back- 
ground configuration changes. 

As was the case in the constant- S physical acous- 
tics analysis, it is convenient to prc-multiply and post- 
multiply the Fresnel equation by S +1 / 2 , thereby re- 
writing the Fresnel equation as 

det {uj 2 I - 2u {v ■ k) I 

- C 2 k 2 - (v ■ k) 2 ij - M 2 } = , (74) 

where again 

C 2 =E 1 ' 2 DE 1 / 2 ; M 2 = -~>l 2 AS 1 / 2 . (75) 

Thus by adopting the eikonal approximation implicit in 
the Fresnel equation one has reduced the number of ma- 
trices one needs to deal with to C 2 and M 2 , now without 
needing to strictly enforce the spatial and temporal con- 
stancy of 5 (or D). The Fresnel equation now becomes 



det 



wo 



2 

k I I 



2;„2 



C z k 



0. 



(76) 



Imposing mode decoupling, which is now equivalent 
to the Fresnel equation factorizing into two quadratics, 
forces M 2 and C 2 to be simultaneously diagonalizable. 
[In which case we recover equation (|46|) .] Finally, impos- 
ing mono-metricity enforces C 2 — c 2 I, with M 2 sym- 
metric and singular but otherwise unconstrained. 

Thus an analysis in terms of the Fresnel equation leads 
to the same conclusions as the direct physical acoustics 
analysis in the special case of constant S, or the special 
case of constant D — these are the same algebraic con- 
straints as were obtained when analyzing equations i|54|) 
and JUnj), albeit in a slightly different physical regime. 



an example of two interacting phonon modes whose wave 
equation is governed by the second-order coupled system 
of PDEs (|T%|) or equivalently (|T§|l , and whose dispersion 
relation is governed by the quartic Fresnel equation lf7T))l. 

It is important to note that fully decoupling the wave 
equation in a completely general background is alge- 
braically more restrictive than the problem of decou- 
pling the Fresnel equation in the geometrical acoustics 
limit. Indeed in the geometrical acoustics limit decou- 
pling places algebraic constraints on the coefficients of the 
wave equation that are equivalent to physical acoustics 
in the special case where the matrix S or the matrix D is 
both position and time independent. This situation is to 
a good approximation relevant at and near the centre of 
the BEC cloud. The subtleties involved in implementing 
decoupling into independent modes is surprisingly more 
complex than one might at first envisage. 

If for the sake of discussion we insert specific numbers 
relevant to a BEC mixture based on hyperfine states in 
87 Rb, we find uiq « 150 kHz. For a harmonic magnetic 
trap this should be compared with a typical trap oscil- 
lation frequency of 100 kHz, though note that for non- 
harmonic traps one can in principle make the trap os- 
cillation frequency arbitrarily small. Similarly we find 
A w 10 29 J, corresponding to a laser temperature of 800 
nK (to be compared to a BEC condensation temperature 
of some fiK). In short, the effects we have been consid- 
ering in this article lie at the cutting edge of present day 
experimental technique. In this regard it is important 
to realise that multi-component BECs have already been 
constructed in the laboratory [3j| • The technological is- 
sues in actually implementing this approach for generat- 
ing the massive Klein-Gordon equation amount to keep- 
ing the background condensates in phase while decou- 
pling the phonon modes in a simple manner. Doing this 
may be experimentally challenging, but there appears to 
be no obstacle in principle to actually implementing the 
model. 

From the general relativity perspective, this article set- 
tles an important matter of principle: It provides an ex- 
ample of an analogue system that can be used to mimic 
a minimally coupled massive scalar field embedded in a 
curved spacetime. Quantum fields of this type are essen- 
tial for any realistic application of analogue spacetime 
ideas to particle physics, and in particular are essential 
for developing condensed matter simulations of quantum 
gravity phenomenology. 



XI. DISCUSSION 

In conclusion, the calculation presented in this article 
is of interest to two separate communities. For the BEC 
community, it provides a specific example of how to tune 
an interacting two-BEC condensate in such a manner as 
to obtain a massive phonon. With the background config- 
urations in phase, but without the fine tuning, it provides 
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APPENDIX A: WAVE EQUATION FOR 
ARBITRARY INITIAL PHASES 



If we start our calculations with a more complex phase 
relationship between the two condensates then the wave 
equation l(TB|) gains additional terms. 

The result we obtain for two arbitrary initial — possi- 
bly even time-dependent — background phases is: 



d t (Z- 1 0) = -d t 1 V Yf> ) - V ( V E 
+ V ( (d-VZ^v) V6M A H 



^density ^ Cphase 



v ( VET 1 a 



phase J 



~l~ Cphasc ^density ^ } 



phase 



••)} 



(Al) 



asc "h ^density ^ 

In addition to those matrices defined in the previous cal- 
culations, three new matrices show up. 

The matrix V simply contains the two background ve- 
locities of each condensate, 



V = 



VAO 
VBO 



(A2) 



now with two possibly distinct background velocities, 

VfiUo, 

(A3) 



VAO = — V0AO, 

m 



VBO 



Additionally, we also obtain two completely new ma- 
trices, which vanish in the case of identical initial phases 



A AB ■= Ao - 0bo = 0. (A4) 
These new matrices are the coupling-phase matrix, 



a 



phase 



A sin Aab 

h 




(A5) 



and the coupling-density matrix 



X sin Aab 



density 




(A6) 



In the case of identical phases equation IjAlfl simplifies to 
the wave equation i|18|) found in the previous calculation. 
Specifically the first two lines of (|A1(I are exactly the same 
as (I18I1 . while the next two lines represent a mass-shift, 
and the last two lines correspond to a damping-term and 
a smoothing term respectively. 
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